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We investigate excitation energies of hyperbolically deformed S = 1 spin chains, which are 
specified by the local energy scale = cosh j A, where j is the lattice index and A is the 
deformation parameter. The elementary excitation is well described by a quasiparticle hopping 
model, which is also expressed in the form of hyperbolic deformation. It is possible to estimate 
the excitation gap A in the uniform limit A — >■ 0, by means of a finite size scaling with respect 
to the system size A'^ and the deformation parameter A. 
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1. Introduction 

A role of numerical study in condensed matter physics 
is to analyze ground state properties and also low-energy 
excitations in the thermodynamic limit. Precise determi- 
nation of the excitation gap is important, since presence 
of a finite excitation gap corresponds to the finitely cor- 
related property of the system. In numerical analyses of 
excitation gap, the finite size scaling-'^'^^ has been em- 
ployed. This is because computational resources are lim- 
ited, and therefore the direct treatment of infinite system 
is occasionally difficult. "^^ 

Density matrix rcnormalization group (DMRG) 
method is widely applied to one-dimensional (ID) quan- 
tum systems, since the method provides precise low- 
lying eigenvalues of large scale systems, up to hundreds or 
even thousands of sites. Although it is possible to apply 
DMRG method to systems with periodic boundary con- 
ditions (PBC),^'^'^^ a majority of applications are per- 
formed under the open boundary conditions (OBC) . This 
is because numerical implementation is much easier with 
the use of OBC. However, presence of boundary correc- 
tions often makes precise scaling analyses difficult. To 
suppress such a corrections, so called the smooth bound- 
ary condition has been applied. It is reported that 
an adiabatic decay of interaction strength toward the sys- 
tem boundary drastically decreases the boundary effect 
on one-point functions, such as the energy density. 

The smooth boundary condition is not applicable when 
one is interested in elementary excitations. This is be- 
cause presence of a small energy scale near the bound- 
ary induces a fictitious excitation localized around the 
boundary. We therefore considered to the opposite direc- 
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tion, and proposed the hyperbolic deformation 
introduced position dependence to ID lattice Hamilto- 
nians, where the energy scale increases toward the both 
ends of the system. As a consequence of this non-uniform 
deformation, the excited quasiparticle is weakly bounded 
around the center of the system. The boundary effect for 
the excitation energy is reduced, since the quasiparticle 



does not reach the boundary. 

In the previous study, we chiefly applied the hyper- 
bolic deformation to the free Fermionic lattice model, 
and showed that a two-parameter scaling function ex- 
ists with respect to the system size N and the deforma- 
tion parameter A.^^^ It is expected that a wide class of 
ID systems under the hyperbolic deformation obeys this 
kind of two-parameter FSS, which could precisely esti- 
mate the excitation gap. In this article we investigate the 
efficiency of this FSS in the context of hyperbolic defor- 
mation, when it is applied to correlated systems, such as 
the = 1 AKLT chain^^^ and the S = 1 Heisenberg chain 
with uniaxial anisotropics (XXZ-I-D). In the next section 
we introduce hyperbolically deformed Hamiltonians for 
these systems. In §3 the form of the two- parameter FSS 
is shortly reviewed. Numerical results on the spin chains 
arc shown in §4. Conclusions are summarized in the last 
section. 

2. Hyperbolic Deformation 

Let us consider a group of ID quantum spin chains, 
which are characterized by the non-uniform Hamiltonian 



H{\) 



L/2-] 

E 

j = -L/2+l 



cosh j A h 



i,i+i 



(2.1) 



where j represents the lattice index, L the system size, 
and A a nonnegative deformation parameter. The oper- 
ator j+i specifies the interaction between neighboring 
sites. An example is the bilinear-biquadratic interaction 

Vi + -^sin6i(Sj -Sj+i)^ (2.2) 



hj .jj_i = JcosO Sj 



between neighboring S = I spins, where J > is the cou- 
pling constant, and where the angle 6 determines the ra- 
tio between bilinear and biquadratic interactions. When 
A > 0, the relative bond strength = cosh j A is posi- 
tion dependent, and the weakest bond is located at the 
center of the system j = 0. The Hamiltonian H{X) when 
A > can be regarded as a deformation — the hyperbolic 
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deformation from the uniform one H{X — 0). 

The case tan6» = 1/3 and J = l/cos6» in Eq. (2.2) 
corresponds to the AKLT interaction.^^' In this case, 
the deformed Hamiltonian is explicitly written as 



H 



AKLT 



(A) = 

L/2-1 

cosh j A 

i=-L/2+l 



(2.3) 



which can be called as the deformed AKLT chain. De- 
spite of the position dependence in _ff^^'^'^(A), the mech- 
anism of the Z2 X Z2 symmetry breaking^^^ is preserved, 
since the Hamiltonian can be written as a linear com- 
bination of local projections. The corresponding ground 
state is the uniform valence-bond-solid state. 

A similar uniformity in the ground state is observed 
for the hypcrbolically deformed 5 = 1 antiferromagnetic 
Heisenberg spin chain, whose Hamiltonian 



H 



AFH 



(A) 



L/2-1 

cosh jA Sj ■ Sj_|_i 

=-L/2+l 



(2.4) 



corresponds to the choice with J — 1 and 6 — in 
Eq. (2.2). Even in the case A > the expectation value 
{Sj ■ Sj_|_i) calculated for the ground state is nearly posi- 
tion independent.-'^*^ Also the deformed transverse- field- 
Ising model, which is defined by the Hamiltonian 

L/2-1 

H^^\X) = Yl co^li-?'^ [^'J + 444+1 + 1^4+1] ' 

j=-L/2-|-l 

(2.5) 

where = (Sj,sJ,s|) represents the 5 = 1/2 spin at 
the jth site, exhibits a similar uniformity deep inside 
the system, even at the critical point F — l}^^ These 
uniform property can be qualitatively explained by the 
path-integral representation of imaginary time evolution 

U = exp —tH{X) under the hyperbolic geometry. 

3. Scaling Form for the Excitation Energy 

Let us consider the infinitely long and uniform 5=1 
Heisenberg spin chain, which corresponds to the double 
limit A and L ^ 00 of H^'^^{X) in Eq. (2.4), as a 
reference system. The model has finite excitation energy 
A from the singlet ground state, which is known as the 
Haldanc gap.^°'^^^ It was reported that the quasiparticle 
picture well holds for the magnetic excitation. When 
only a quasiparticle is present, the energy dispersion of 
the quasiparticle can be described by the effective Hamil- 
tonian 

^eff ^ _^ ^ (gt g .^^ + -t^^g + (2t + A) ^ c]cj 

3 3 

(3.1) 

in the vicinity of zero momentum, where t is the hopping 
parameter, and where cj and c^, respectively, represent 
creation and annihilation of the quasiparticle. The eigen- 
value of the zero-quasiparticle vacuum is trivially zero, 
and that of the one-quasiparticle state is larger than zero 
when A is positive. If more than two quasiparticles are 



present, interaction terms should be included into the 
effective Hamiltonian.^'^' ^''^ 

We conjecture that such a quasiparticle picture 
also holds for the hypcrbolically deformed Hamiltonian 
J^^^^(A) in Eq. (2.4), and that the effective Hamiltonian 
can be written as 



eff 



L/2-1 

(A) = coshjA 

i=-L/2+l 



(4' 



+ (2t + A) (c] 



Cj + c\ 



■]+i^j+i 



(3.2) 



in the form of hyperbolic deformation. 

As a preparation for checking the correspondence be- 
tween if^^^(A) and H^^{X), the authors investigated 
one-particle state of in our previous study. ^"^^ 

The lowest-energy one-particle state is a shallow bound 
state, where the quasiparticle does not reaches to the 
both end of the system when A is sufficiently large. The 
lowest eigenvalue, which we denote as Aj;^(A), converges 
to A in the limit L ^ 00 and A ^ 0. We found that the 
correction t^i^(X) — A satisfies the two-parameter scaling 



2 A^(A) 



A 



t 



= 9 



{L + lfX 



(3.3) 



with respect to the system size L and the deformation 
parameter A. The scaling function g[y\ satisfies (?[0] = tt^ 
and g{y ^ 1] ~ y/\/2. We conjecture that this scal- 
ing form is also applicable for elementary excitation of 
gapped spin chains in general. 

4. Two-parameter Scaling Analysis 

We perform numerical analysis of the hyperbolically 
deformed spin chains defined in §2, by means of the 
DMRG method. In order to avoid the quasi four-fold 
degeneracy,^'') we put 5 = 1/2 spins at the both ends 
of the system; these boundary spins are not counted 
when wc refer to the system size L. For the bond between 
boundary 5=1/2 spin and the neighboring 5 = 1 spin, 
we set weaker interaction strength J^^^ = 0.5088. Actu- 
ally the value of J^,^^ is not relevant to the elementary 
excitation and its energy if XL is sufficiently large. We 
keep at most m = 100 block spin states in the DMRG cal- 
culations. Numerical convergence in finite system sweep- 
ing is accelerated by use of the wave function prediction 
method. ^^~^*) We explain numerical details in the ap- 
pendix A. 

Concerning to the ground state of the hyperbolically 
deformed 5=1 spin chains that we have introduced, 

the z-component of the total spin is zero. Let us denote 
the corresponding ground state energy by E^f \X). We 

also calculate the lowest eigenvalue E^j^\x) when the z- 
component of the total spin is one. The magnetic excita- 
tion energy is then expressed as their difference 

A^{X)=E'i\x)-Ef\x). (4.1) 

We have used the same notation A]^{X) that appears in 
Eq. (3.3), since we expect that A^{X) in Eq. (4.1) also 
satisfies the scaling relation in Eq. (3.3). 
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Fig. 1. Scaling plot for the hypcrbolically deformed 5=1 Heiscn- 
berg chain. Best fit is realized when A = 0.410485 and t = 7.381. 



Fig. 3. Scaling plot for the hypcrbolically deformed XXZ chain, 
where a = 0.8 and D = in Eq. (4.2). The parametrization 
A = 0.20502 and t = 13.01 is used. 
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Fig. 2. Scaling plot for the hyperbolically deformed AKLT chain. 
The parametrization A = 0.7002483 and t = 0.51542 is used. 
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Fig. 4. Scaling plot for the case a = 1 and D = 0.2 in Eq. (4.2), 
under the parametrization A = 0.288240 and t = 10.77. 



For H^^^ {X) in Eqs. (2.4), we confirm the presence 
of scaling function g under the choice A = 0.410485 and 
t = 7.381 after some trials of determining these param- 
eters. Figure 1 shows the scaling result for A^(A) — A 
when L = 100, 200, and 500. These data agrees with 
A^(A) — A shown by solid curve, which is drawn from 
the effective one-particle model H'^^{X) in Eq. (3.2). We 
also calculate ^i^{X) — A for the deformed AKLT chain 
defined by ^^^^'^(A) in Eq. (2.3), and show the result 
in Fig. 2. For this case the best fit to H'^^{X) is real- 
ized when A = 0.7002483 and t = 0.51542. For both 
Heisenberg and AKLT chains, the quasiparticle picture 
well holds even under the hyperbolic deformation, and 
the estimated gaps A by these finite size scalings are 
consistent with known values.^^'^^'^^'^*^^ 

Let us check the validity of two-parameter scaling on 
another model, the S = 1 XXZ spin chain with uniaxial 
anisotropy. Under the hyperbolic deformation, the corre- 
sponding Hamiltonian is written as 



H 



XXZ 



(A) 



L/2-1 

cosh j A 

=-L/2+l 



qX qX I qY qY 



L/2 

+ D cosh(i-i)A [S] 



(4.2) 



Figure 3 shows the scaling plot when a = 0.8 and D = 0. 
For this case the scaling parameters are determined as 
A = 0.20502 and t = 13.01. Figure 4 is the case when a = 
1 and D = 0.2. For this case we obtain A = 0.288240 and 
t = 10.77. Again the two-parameter scaling well holds 
for H-^^^{\), and the obtained value of A agrees with 
previous studies. ^^"^^^ 

A profit of the two-parameter scaling is that the ef- 
fective hopping parameter t is obtained simultaneously 
with the gap A. It is well known that the spin velocity 

V = V2tK (4.3) 

is essential for the thermodynamic property of the spin 
chains at low temperature. (See Appendix B.) From 
the data in Fig. 1, the velocity at the isotropic Heisen- 
berg point is obtained a,s v = 2.462, which is consistent 
with V = 2.4691 reported by S. Todo.^^) At the AKLT 
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point, the estimated value of the velocity from Fig. 2 is 
V = 0.84962. The value deviates from v = 0.825 reported 
by K. Okunishi et al,^^^ and the clarification of this dif- 
ference is a remaining problem. 

5. Conclusion 

We have investigated the elementary excitation en- 
ergy of the hyperbolically deformed S = 1 spin chains. 
Two-parameter scaling for the excitation energy A^(A) 
coincides with the effective one-q\iasiparticle picture in 
Eq. (3.2), which is also written in the form of hyperbolic 
deformation. The fact makes it possible to estimate both 
the bulk gap A and the effective hopping parameter t ac- 
curately by introduction of the hyperbolic deformation. 

There are some variations of the non-uniform defor- 
mation applied to spin chains. A simple example is the 
exponential deformation introduced by Okunishi. ^^■'^^^ 
Similar to the hyperbolic deformation, the correlation 
length becomes finite when the site-dependent energy 
scale fj = e^^ is introduced. The behavior of the excited 
quasiparticle from the gapped ground state under this 
exponential deformation is our interest for future study. 

Acknowledgements 

The authors thank to Koichi Okunishi for valu- 
able discussions. This work was partly supported by 
Grants-in-Aid for JSPS Fellows, Grants-in-Aid from 
the Ministry of Education, Culture, Sports, Science 
and Technology (MEXT) (No. 19540403, No. 20340096, 
and No. 22014012), and the Global COE Program 
(Core Research and Engineering of Advanced Materials- 
Interdisciplinary Education Center for Materials Sci- 
ence), MEXT, Japan. 

Appendix A: Numerical Details 

We explain several numerical details on the calcula- 
tion of £;["'(A) and £;[^^(A) in Eq. (4.1) by means of 
the DMRG method. To simplify the formulation without 
loosing generality, we consider H{X) defined by Eq. (2.1). 
According to the custom in DMRG, let us split H{X) into 
three parts 



HiX) = H^iX)+ho.+H^{X)., 



(A-1) 



where Hj^{X) and J?jj(A) are, respectively, the left and 
the right block Hamiltonians, which are defined as 

-1 

HlW = coshjA hjj^i 

j=-L/2+l 
L/2-1 

Hr{X) = Yl coshjA /I,. ,.+1 . (A.2) 

These block Hamiltonians satisfy a recursion relation^^^ 
= cosh A Q — /i_2 _i 
+ 2 cosh X HI- HI* {X) 

— cosh A hi 2 — ^2,3 

+ 2 cosh XH%- H*j^ (A) , (A-3) 



where H£ and H'^ represent the block Hamiltonians for 
the {L — 2)-site system 



Hl{X) 



Y cosh(j 

-L/2+1 



1)A 



L/2-1 

H*r{X) = Y co^l^O' - kj+i ' (A-4) 

and where H^* and H^ represent those for {L — 4)-site 
systems 

-2 

^r(A) = Yl oosh{j-2)Xh^j+i 

j=-L/2+l 



Hr{X) 



L/2-] 

E 

j=2 



cosh(j - 2)A hjj+i . (A-5) 



The recursion relation Eq. (A. 3) is useful when one per- 
forms the infinite system DMRG iteration for the prepa- 
ration of the block Hamiltonians to start the following 
finite system DMRG sweeps. It is also possible to per- 
form the infinite system DMRG method for a uniform 
spin chain A = 1, and increase A afterword adiabatically 
during the following finite system DMRG sweeps. 

A special care, which should be taken for the hyperbol- 
ically deformed system, is the subtraction of the energy 
expectation value from the Hamiltonian. Since the factor 
cosh J A increases rapidly with respect to the value of 

E^f ^ becomes huge when -^A is large. This is a cause of 
numerical error. Thus we shift the origin of the energy 
so that the ground-state energy becomes almost zero. 
This energy shift can be performed during the finite size 
sweeping. First we obtain the ground state j^'g) diago- 
nalizing the Hamiltonian H{X), which is represented as 
the super-block Hamiltonian in DMRG formulation. We 
calculate the expectation value {hj jj^i) for each bond 
between the active two sites between the left and the 
right blocks. Every time we calculate {hjj_^i), we sub- 
tract coshXj (hj j_^_i) from the bond Hamiltonian, and 
create the renormalized block Hamiltonians after this 
subtraction. Within several sweeps the numerical error 
in {hj j_^_i) almost vanishes, and the ground state energy 
calculated for the Hamiltonian under the energy shift 



HiX) 



L/2-1 

{H{X)) = Y '^os^-?'^ 

j=-L/2+l 



(A-6) 

becomes zero within the tiny numerical error. When we 
put S = boundary spins, we also perform this energy 
shift for the additional boundary interactions. 

Appendix B: Spin Velocity 

We shortly review a physical interpretation of the spin 
velocity v = \/2tK in Eq. (4.3). In special relativity the 
energy E of a particle that has momentum p and mass 
m is given hy E ^ \J p^c? + m^c^, where c is the velocity 
of light. Expanding by p one obtains the non relativistic 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 5 



approximation 

E mc^ +p^/2m = mc^ +^^0^/27X10^. (B-1) 

On the other hand, the dispersion of H'^^ in Eq. (3.1) is 
given by E^^ = A + 2t(l — cos A;), where k is the wave 
number of the quasiparticle. Expanding by k one obtains 

(B-2) 



£;eff ^ A + ifc^ = A + fc^2iA/2A. 



It should be noted that both the lattice constant and 

k arc dimcnsionlcss. Thus mc^ corresponds to A, and p 
corresponds to k. Comparing Eq. (B-1) and Eq. (B-2), 
one finds that the light velocity c corresponds to a di- 
mcnsionlcss velocity V ~ ^/2tA, which is nothing but the 
velocity of the spin wave when the system is gapless. 
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